Quantum gates and multi-particle entanglement by Rydberg excitation blockade and 

adiabatic passage 
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We propose to apply stimulated adiabatic passage to transfer atoms from their ground state 
into Rydberg excited states. Atoms a few micrometers apart experience a dipole-dipole interaction 
among Rydberg states that is strong enough to shift the atomic resonance and inhibit excitation of 
more than a single atom. We show that the adiabatic passage in the presence of this interaction 
between two atoms leads to robust creation of maximally entangled states and to two-bit quantum 
gates. For many atoms, the excitation blockade leads to an effective implementation of collective- 
spin and Jaynes-Cummings-like Hamiltonians, and we show that the adiabatic passage can be used 
to generate collective Jx = G eigenstates and Greenberger-Horne-Zeilinger states of tens of atoms. 

PACS numbers: 03.67.Bg, 03.65.Ud, 42.50.Dv 



Entanglement is a central property of quantum me- 
chanical systems and an important resource for appli- 
cations in quantum information and quantum metrology. 
Entanglement is subject of broad theoretical interest and 
experimental implementations in various systems have 
been demonstrated. One scheme for creation of entan- 
glement and quantum gates between neutral atoms uti- 
lizes the large dipole moments of highly excited Rydberg 
states resulting in a dipole-dipole interaction which is 
strong enough to shift the atomic energy levels and pre- 
vent more than one atom from being excited to the Ryd- 
berg state by a resonant laser field [l|, ■ In alkali atoms 
Rydberg states with principal quantum number ri > 70 
have lifetimes > 100 /iS and experience a dipole-dipole in- 
teraction strength, E/h, above 100 • 27r MHz when atoms 
are separated less than 5 /xm The accompanying 
suppression of excitation for n up to 80 [1, Q and the in- 
fluence of the excitation blockade on coherent collective 
dynamics 0] was observed in cold gases. 

We propose to combine the Rydberg blockade mecha- 
nism with the rapid adiabatic laser pulse sequence known 
from Stimulated Raman Adiabatic Passage (STIRAP) 
0, H, [l3|- We consider atoms with two lower levels 
and a Rydberg level in a ladder structure coupled by 
two resonant laser fields with Rabi frequencies fJi and 
Vlr as shown in Fig. [TJa). The STIRAP process ap- 
plies the "counter intuitive" pulse sequence of Fig. [H^b) 
to transfer population from the ground state |1) to a 
highly excited Rydberg state \r) by adiabatically follow- 
ing ll| a dark state which never populates the inter- 
mediate state 1 2). If we assume a relative phase, (j)r{t), 
between the Rabi frequencies f2i and f2r, a single atom 
exposed to the fields follows the dark state superposition, 



\Di) = cos0|l) - sin^ 



where tan 9 = fli/flr ex- 



presses the relative strength of the two laser fields. Since 
the dark state has zero energy, there is no dynamic phase 
accumulated during the process, but if (pr varies, |Ui ) ac- 
quires the geometric phase, 7i = ~ / sin^ 9d(j)r 12 



controlled geometric phase shift on state |1) can be imple- 
mented if the phase difference between the fields evolves 
between the STIRAP pulses and a second inverted set 
of pulses, shown in Fig. [Ijc). For atoms with a further 
logical state |0) such a pulse sequence enables robust ge- 
ometric phase gates and other one-bit gates [l^l ■ 
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FIG. 1: (a) Three-level ladder system and two laser fields with 
Rabi frequencies Qi and Jlr. (b) STIRAP pulse sequence. 
The FWHM of each pulse is a and the delay between pulses 
within one process is At. (c) Pulse sequence consisting of two 
STIRAP processes separated by AT in time. 

With more atoms the Rydberg blockade comes into 
play, but we may still identify dark states. For two atoms, 
initially in the product state |11) and subject to the same 
interaction with the two resonant laser fields fli and 
the evolution preserves the symmetry under interchange 
of atoms, and hence it is sufficient to consider the Hamil- 
tonian in the symmetric two-atomic basis {|ll),^(|lr)-|~ 

|rl)). J.(|12) + |21)). |r,-), A-fflr) + \r2)). |22)}, 
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where E denotes the energy shift of the state \rr) due to 
the dipole-dipole interaction. This Hamiltonian has one 
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dark state for the two-atom system 

\D2)= . i^ =[(cos^0-sm^g)|ll) (2) 

\/cos4 6l + 2sin^6' 

-cos6lsin6'e*'^'-(|lr) + |rl)) + sin^ 6*122)] . 

We first assume e*'^'' = 1 and apply the counterintuitive 
pulse sequence of Fig. [TJb) . Initially cos = 1 and the 
system is in \D2) = Adiabaticity ensures that we 

remain in \D2) and after the pulses sinO ~ 1 and the 
system is in 



1 



|11) + |22)). 



(3) 



Equation ^ shows, that while the STIRAP process in 
the single atom case ensures that |2) is never populated, 
due to the Rydberg blockade the pair of atoms are adia- 
batically steered into a state populating |22). Moreover, 
([3]) is a maximally entangled state of the two atoms, gen- 
erated robustly irrespective of the precise pulse shapes, 
field strengths and the precise value of the Rydberg in- 
teraction energy. Figure[2Ua) shows the evolution of pop- 
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FIG. 2: (a) Time evolution of the population in 
^(|lr) + |rl)) and |22). (b) Fidelity of the entanglement 
creation as a function of pulse width for different values 
of dipole-dipole interaction, E. \tpf) was found by propa- 
gating the Schrodinger equation ([T]) with initial state 
We show the norm square of the overlap with the target 
state \D2)- The laser pulses are modeled by sin^-pulses, 

0,j{t) = rJmaxj- sin^( ^^*^J°^^ ) for tsj <t< taj + 2(T, where tsj 
is the instant of time when the pulse starts and a the FWHM. 
The simulations are made with parameters ilmax,i/'2iT = 
^max,2/'2TT = 10 MHz, At = 1.1 jis, lifetime of the \r) state, 
Tr = ioO ^ts and in (a) a = 1.5 fis and E/h = 100 • 2tv MHz. 

ulations for realistic experimental parameters, obtained 
from a numerical propagation of the Schrodinger equa- 
tion (see caption) . The decay due to the finite lifetime of 
the |r) states, populated during the process, is incorpo- 
rated as a decay out of the system, whereas the |1) and 
1 2) states are treated as long lived hyperfine sublevels. 
This implies that if the levels are coupled by single pho- 
ton transitions, the field coupling |1) and |2) will have 
a frequency in the radiofrcquency range, while the laser 
field coupling |2) and |r) has a wavelength w 300 nm for 



Cs and Kb. Alternatively the couplings can be obtained 
with two-photon transitions involving four optical fields. 
By appropriately choosing the laser detunings to com- 
pensate for the time-dependent AC Stark shifts during 
the pulses, the corresponding five level system of equa- 
tions with two intermediate optically excited states can 
be effectively reduced to the present three level ladder 
system. 

To investigate the criterion of adiabaticity and the role 
of the Rydberg interaction energy E, in Fig. [Hb), we 
show the fidelity of the creation of the entangled state, 
F ~ |(D2|'0/)P, where \ipf) is the final state calcu- 
lated by propagation of the state vector with the time- 
dependent Hamiltonian ([1]). The simulations show that, 
as long as E is sufficiently large to block the population 
of states with more than one Rydberg excitation, the ex- 
act value of E is not important. With Rabi frequencies 
f^max,i/27r ^ ri,nax,r/27r = 10 MHz, a Rydberg energy 
shift oi E/h= 50 • 27r MHz is sufficient. The time window 
where \r) is populated is determined by the pulse width 
and it is desirable to use the smallest possible width that 
does not violate adiabaticity, yielding a total time of en- 
tanglement generation of 3-4 fis, which is short compared 
with the radiative lifetime of the highly excited Rydberg 
state > 100 ^s. 

The two-atom entanglement scheme can be modified 
to create a two-qubit phase gate. To this end we apply 
a second STIRAP sequence with phase shifted Rabi fre- 
quencies and the pulses in reversed order, cf. Fig. [Ijc), 
that transfers the population back to the |11). With a 
non- vanishing relative phase (t>rit) between fli and fir 
the dark state ^ acquires the geometric phase 
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When wc supplement the atomic level scheme with an- 
other qubit state |0), which is uncoupled from the STI- 
RAP pulses, the gate performed by the two STIRAP pro- 
cesses amounts to multiplication of all two-bit register 
states by phase factors corresponding to a controlled two 
qubit phase gate with phase Aip = 72 — 271. 72 is only 
acquired when the pulses overlap, while 71 is acquired in 
between the two pulse sequences, and A(p can, e.g., be 
controlled by adjusting AT. 

We now show that when more than two atoms in |1) 
are subject to the STIRAP pulse sequence they also be- 
come entangled. Provided all atoms are localized within 
a region of a few /zm, the transition towards the Rydberg 
states is restricted to the coupling of collective states with 
either zero or one Rydberg atom and we implement this 
by truncating the basis so it only includes states with zero 
or one atom in the |r) state and write the symmetric ba- 
sis states of the system as {|ni, n2 = iV — ni, 0), |ni, 712 = 
iV — 7ii — 1, 1)}, where N is the total number of atoms, 
Til and ri2 the number of atoms in |1) and |2) and 
or 1 indicates if \r) is populated with zero or a single 
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atom. The interaction with the radiation fields fli and 

^r, J2f=i-li^ii'tMj{M + ^rit)\r)j{2\+h.c.), can now 
be rewritten 



H{t)^HjM+Hjc{t), 



(4) 



with variable strengths, representing the coupling by the 
fields driving the lower and the upper transition, respec- 
tively. The dynamics of the lower levels can be rewritten 
in a collective spin description, and in the accompanying 
Schwingcr oscillator description, 

Hj^t) = -hVLi{t)J^{i) = -]^hni{t){a\a2 + ai4), (5) 

where a^^^ are annihilation (creation) operators, with 
conventional oscillator commutator relations, for the 
number of atoms in The upper transition is described 

by 



Hjc{t) 



1 



ala 




(6) 



t H 



FIG. 3: (color online) Eigenvalues of H{t)/h of (gl for iV = 6. 
Black solid curves show the results of a numerical diagonaliza- 
tion of the full Hamiltonian. The squares on the left show the 
analytical eigenvalues of Hjc{t)/h and the bullets and trian- 
gles on the right show the analytical eigenvalues of Hj^{t)/h 
for zero and one Rydberg excitation, respectively. Parame- 
ters used are f2maa;,i/27r = 10 MHz, ^Imax.r /2n = 10 MHz, 
a = 1.5 /is. At — 1.1 /is 



where we use Pauli matrices, cr"*" and a~ , to describe the 
transfer of population between states with zero and one 
Rydberg excitation in the atomic ensemble. This Hamil- 
tonian is the quantum optical Jaynes-Cummings (JC) 
Hamiltonian, introduced originally to describe the reso- 
nant interaction between a two-level atom and quantized 
light [14|. JC dynamics has been implemented in strong 
coupling cavity QED experiments, where non-classical 
states of light, such as Fock states and quan tum super- 
position states, are produced efficiently [l5|. Hjc also 
describes the motion of laser driven trapped ions where 
it has been used to generate various non-classical states 
3 ill- Following these proposals, JC dynamics is suf- 
ficient to produce a variety of interesting states of an 
atomic ensemble by the Rydberg blockade. Because of 
its robustness, however, we shall proceed here, and ex- 
plore the states resulting from the STIRAP protocol. 

The Hamiltonian in ([4]) has a dark state with zero val- 
ued energy eigenvalue throughout the STIRAP process. 
This follows because the "parity" operator Q = (—1)"^, 
which inverts the sign of the operators 02 and Oj, anti- 
commutes with H, HQ = —QH. Any cigenstatc of 
H with energy eigenvalue £, then has a partner, &\ip) 
with eigenvalue -£ and the energy eigenvalue spectrum 
is symmetric around zero. The number of eigenstates is 
odd because N atoms induce TV -I- 1 different states with 
no atoms in \r) and N states with one atom in |r), and 
hence there must always be a state with eigenvalue zero. 
The full curves in Fig. [3] show the 13 energy eigenval- 
ues for 6 atoms found from a numerical diagonalization 
of Eq. (lU, and they clearly confirm the existence of the 
dark state throughout the pulse sequence. 

The STIRAP process starts with Oi = and hence 
initially H{t) = Hjc{t) = -~^hnr{t){a2(T+ + a^a') with 
the eigenvalues {0, ±lhD,r{t)y/n2 + 1}, as shown with the 



squares in Fig.[31 The dark state is \ni = N,n2 = 0, = 
0) with all atoms in state |1). Adiabaticity ensures that 
we remain in the dark state of the full Hamiltonian and 
when f2i is turned on and fir reduced, the system ends 
up in the dark state of H{t) = Hj^{t) = ~Mli{t)Jx- 
In this state, the system can have either no or a single 
Rydberg excitation leaving K = NotK~N — \ atoms 
in the |1) and |2) states. The eigenvalues of —Q,i{t)Jx are 
-VLi{t){-K/2, + 1, K/2 - 1, K/2} and £ = 

occurs for K even. Figure [3] shows the eigenvalues of 
—^i{t)Jx when \r) is not populated (bullets) and when 
one atom is excited to \r) (triangles) for = 6. For 
even, the dark state does not populate |r), while for A^ 
odd, the final dark state is the state with one Rydberg 
excitation and A" — 1 atoms in the J^, = eigenstate. In 
general we write the final dark state 



Dn) = 



0) 

= 0) ® \r))sy^ 



if A^ is even 
if A^ is odd 



(7) 



where {■)sy7n indicates that the state is symmetrized with 
respect to |r), such that any atom is Rydberg excited 
with equal weight. 

The STIRAP protocol produces a. \Jx = 0) multi- 
particle entangled state, and precisely this state reaches 
the Heisenberg limit of phase sensitivity in entanglement 
enabled precision metrology [l^, [l^l- If the two lower 
states are, e.g., the hyperfine states of the Cs clock tran- 
sition, the presented entanglement scheme thus consti- 
tutes an ideal preparation of the system for an atomic 
clock. 

The state ([7]) has none or a single Rydberg excited 
atom, depending on the number of atoms initially in the 
|1) state, ni, being even or odd. Following [2^, this can 
be used to prepare a Grcenberger-Horne-Zcilinger (GHZ) 
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state 2l| of the system. We first prepare all our atoms 
in a spin coherent state, i.e., a product of superpositions 
of two ground states |0) and 




na,ni) 



(8) 



where ni is the number of atoms in the state |1) and 
no ^ N — rii. Applying the pair of STIRAP processes in 
Fig. [He) to the entire system, transfers each component 
|no,ni) to a state with none or a single Rydberg exci- 
tation and back. An energy shift of the Rydberg state 
or a phase shift of the laser coupling the Rydberg state 
can now be used to provide states populating the Ryd- 
berg state with a phase, i = e^^^^. We can write this 
phase, e*''^"!) = (e*'^/'* -I- (-l)"ie-*'^/^)/\/2 leading after 
the second STIRAP process to the final GHZ state. 



N 



ni=0 

/|0) 




e^''("i)|no,7ii) (9) 
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When N is increased, to maintain adiabaticity one 
must require longer pulses and/or a stronger Rydberg 
blockade. As indicated by the energy spectra in Fig. [31 
there are crossing regions, where special care should be 
taken, and we anticipate that control theory may be 
used to find optimal pulse shapes. Using the simple 
sin^-pulses, we have investigated the preparation of the 
\Jx =0) and GHZ ^ states for up to TV = 10 tak- 
ing into account the coupling to states with two Ryd- 
berg excited atoms and the accompanying energy shift. 
Note that the explicit Rydbcrg-Rydberg interaction en- 
ergy does not anti-commute with (— 1)"=^ and the spec- 
trum is not symmetric around zero to this order of pre- 
cision - the non-resonant coupling to the states with the 
large interaction energy shift causes a perturbation of 
the energy levels in the figure. This perturbation im- 
plies a small dynamic phase on the dark state, even if 
adiabaticity is not violated. We find a population of the 
\Jx =0) or GHZ state above 0.995 for peak Rabi frequen- 
cies of 10 • 27r MHz, a Rydberg interaction, E/2Tr = 400 
MHz and pulse widths of 50 /xs. These numbers require 
the use of Rydberg levels with n > 100 where a lifetime 
above 1 ms can be achieved for atoms in a cryogenic envi- 
ronment . If higher Rabi frequencies can be achieved, 
adiabaticity can be maintained with shorter pulses and 
hence Rydberg states with shorter lifetimes can be used. 

In conclusion, we have demonstrated that the Rydberg 
excitation blockade mechanism in conjunction with rapid 
adiabatic passage processes provide rich opportunities to 
prepare two- and multi-atom entangled states with con- 
fined samples of atoms, e.g., in optical dipole traps or 
small lattice arrays. Our calculations indicate the possi- 
bility of high fidelity generation of entangled states and 



quantum superpositions states with tens of atoms, and 
we propose to apply control theory methods to optimize 
pulse shapes and reach even larger systems. 
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